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Abstract. A surface model of Nambu and Goto is studied statistical mechanically by using the canonical 
Monte Carlo simulation technique on a spherical meshwork. The model is defined by the area energy 
term and a one- dimensional bending energy term in the Hamiltonian. We find that the model has a large 
variety of phases; the spherical phase, the planar phase, the long linear phase, the short linear phase, the 
wormlike phase, and the collapsed phase. Almost all two neighboring phases are separated by discontinuous 
transitions. It is also remarkable that no surface fluctuation can be seen in the surfaces both in the spherical 
phase and in the planar phase. 

PACS. 64.60.-i General studies of phase transitions - 68.60.-p Physical properties of thin films, nonelec- 
tronic - 87.16.Dg Membranes, bilayers, and vesicles 



. 1 Introduction 

' Membranes have long been interested in biological physics 
and in softmatter physics [Tl[51[3] . The shape of membranes 
has been studied experimentally and theoretically 4,5,6, 

, WSS\ ■ The statistical mechanical studies for the conforma- 

■ tion and the elastic property of membranes were reviewed 
;in 0. 

The curvature surface model of Helfrich and Polyakov 
' pUlimiT^ plays a role for describing the shape of mem- 
, branes from the two-dimensional differential geometrical 
• or the classical mechanical view point. In fact, the surface 
[ shapes are understood as the conditions that minimize the 
I curvature energies of the models such as the area difference 
' bilayer model and the minimal model [51171ISI ■ It is also in- 
, teresting to study how thermal fluctuations influence and 

■ change the membrane shape [T31IT4]. The prolate-oblate 
^ transition was in fact observed as a phenomenon driven by 

thermal fluctuations [13]. Therefore, it is still worthwhile 
to study numerically the shape transformation phenom- 
ena as phase transitions in the surface models. 

Contrary to the shape transformation, the surface fluc- 
tuation phenomena were extensively studied theoretically 
[15l[16l[ll] and numerically [l8lfT9l[20l[2T|[22ll23p4l[25] in 
the surface model of Helfrich and Polyakov. It was re- 
ported recently that the transition is of first-order on the 
fixed connectivity spherical surfaces [26ll27l[28] . We see 
that the transition is accompanied by the collapsing tran- 
sition, which can be considered as a shape transformation 
transition between the smooth phase and the collapsed 
phase. 

Here we comment on the terminology used in this pa- 
per for membranes models. The Gaussian bond potential 



is defined by the summation over the bond length squares 
on triangulated surfaces, while the Nambu-Goto energy 
is defined by the summation over the area of triangles. 
Both potentials play a role of the microscopic surface ten- 
sion in the surface models. A surface model defined by 
the Gaussian bond potential with an extrinsic curvature 
energy term is called the Helfrich and Polyakov model, 
while the one defined by the Nambu-Goto energy is called 
the Nambu and Goto model. On the contrary, the fluid 
surface model is the one defined on dynamically triangu- 
lated surfaces, and hence the vertices diffuse freely on such 
fluid surfaces. The fixed-connectivity (or tethered) surface 
model is defined on fixed connectivity triangulated sur- 
faces, and hence the vertices can only fluctuate locally on 
such fixed-connectivity surfaces. 

We must recall that the surface models for describ- 
ing the surface fluctuation phenomena can be classified by 
the bending energy (extrinsic curvature or intrinsic curva- 
ture) , the surface tension energy (Gaussian bond potential 
or Nambu-Goto energy), and the fluidity (flxed connec- 
tivity or fluid). Surface models show a variety of shapes, 
which strongly depends on these physical quantities [121 
|30] . The surface model of Nambu-Goto with an intrinsic 
curvature energy has a variety of phases even on fixed 
connectivity spherical surfaces 

One additional element that characterizes the surface 
model is the cytoskeletal structure [21], where the bend- 
ing energy is not always a two-dimensional one |321l33j . 
The cytoskeletal structure is a lattice analog of cytoskele- 
tons or microtubules in biological membranes. The one- 
dimensional bending energy, which is defined on the skele- 
tons of the surface, can also give the mechanical strength 
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to the surface for maintaining the shape. We can call 
a surface model with such cytoskeletal structure as the 
compartmentalized model. It was shown that such com- 
partmentalized models on fluid surfaces have a variety of 
phases including the planar phase and the tubular phase 

Therefore, it is very interesting to study a surface model 
defined by combining these elements such as the Nambu- 
Goto energy and the cytoskeletal structure. In this pa- 
per, we study the Nambu-Goto surface model with a one- 
dimensional bending energy on fixed connectivity spher- 
ical surfaces. Although the lattice is a triangulated one 
and is identical to that for the standard model without 
the compartmentalized structure, the lattice is considered 
to be the compartmentalized one. In fact, the lattice is ob- 
tained from the finite-size compartmentalized one in the 
limit of L ^ 1 (in the unit of bond length), where L is the 
length of the compartment boundary between two junc- 
tions. The compartment size can also be characterized by 
n the total number of vertices in a compartment. The one- 
dimensional bending energy is defined on the compart- 
ments of the model of finite n [32] ■ In Ref . [35] . we stud- 
ied a compartmentalized model with the one-dimensional 
bending energy in the limit of n — > and reported certain 
differences between the results of the compartmentalized 
model at n ^ and those of the conventional model in 
[27] . It is remarkable that the results in this paper are 
significantly different from those in [35]. Only difference 
between the model in this paper and that in [35| is in the 
bond potential; the Nambu-Goto energy is assumed in the 
model of this paper while the Gaussian bond potential is 
assumed in I35l . 



2 Model and Monte Carlo technique 

Triangulated lattices are obtained by splitting the icosahe- 
dron, which has 12 vertices of coordination number g = 5. 
By dividing every edge of the icosahedron into £ pieces, 
we have a triangulated lattice of size + where 

the generated TV— 12 vertices are of coordination number 

g=6. 

The surface model in this paper is defined by the Hamil- 
tonian S, which is given by a linear combination of the 
Nambu-Goto energy Si and the one-dimensional bending 
energy S2 such that 

S = Si + bS2, = 52 = ^(l-t,.t,), (1) 

A i] 

where b is the bending rigidity. in Si is the area of 
the triangle A, and in S2 is a unit tangential vector of 
the bond i. in Si is the sum over all triangles A, and 
J2ij is the sum over bonds i and j. 

In S2, the combination 1— - tj of the vectors and tj 
is defined as follows: Three combinations are defined at the 
q = 6 vertices, and 5/2 combinations are defined at the q = 
5 vertices. Figures IHa) and Wlh) show the combinations 
1 — tl • t2 and 1 — tl • (t2-l-t3)/2, which are respectively 
typical of the q — 6 vertices and of the q~5 vertices. The 



combination 1— ti • t2 is included in the summation of S2 
with the weight of 1, while the combination 1— ti • (t2-l-t3)/2 
is included in S2 with the weight of 1/2. This is the reason 
why the 17 = 5 vertices have 5/2 combinations. 




(a) (b) 

Fig. 1. (a) A combination 1— ti ■t2 is included in the bending 
energy with the weight of 1 at the g = 6 vertices, and (b) a 
combination 1— ti ■ (t2-|-t3)/2 is included in the bending energy 
with the weight of 1/2 at the g = 5 vertices. 

The partition function of the model is defined by 

Z= \{dX,e^^[-S{X)], (2) 
•' i=i 

where A'i(e R^) is the position of the vertex i on the 
triangulated spherical surface. 

We should note on the problem of well-definedness of 
the Nambu-Goto surface model. Our model in this pa- 
per will be included into a class of well-defined models of 
Nambu-Goto. The surface model with Hamiltonian S — Si 
is known to be ill-defined [36]. In fact, it is easy to see nu- 
merically that the surface is just like a chestnut bur in the 
model defined by S — Si. Since the spines are extremely 
long, the resulting surface appears to be a gas. Some of 
the vertices are spread over large space in R'^. This ill- 
definedness persists even when the Hamiltonian includes 
the standard bending energy of the form 1— rii-rij , where 
is a unit normal vector of the triangle i. On the contrary, 
the model turns to be a well-defined one if certain types of 
bending energies are included in the Hamiltonian |28p29j. 
A bending energy, which is an extrinsic curvature, makes 
the model well-defined, and it was shown that the model 
has the smooth phase and the collapsed phase [28] . It was 
also reported that an intrinsic curvature energy makes the 
model well-defined, and a variety of phases can be seen in 
the model [25j . 

The three-dimensional integrations in the partition func- 
tion are simulated by the random three-dimensional shift 
of X in the canonical Monte Carlo (MC) simulations. The 
new position X' is given by X' — X+5X , where 5X is ran- 
domly chosen in a small sphere. The radius of the small 
sphere is fixed at the beginning of the simulation to main- 
tain about 50% acceptance rate. 

The surface size assumed in the simulations is iV = 
1442, N = 2562, and N = 4842. The total number of 
Monte Carlo sweeps (MCS) after the thermalization MCS 



H.Koibuchi; Shape transformation transitions of a tethered surface model 

(a) (b) 



is about 1x10^ ~ 2x10^ for the surfaces of size iV=1442, 
2x10® - 3x10® for the surfaces of size iV = 2562, TV = 4842 
at b close to the transition point, and relatively smaller 
number of MCS is assumed at b far from the transition 
point. 




(c) 

Fig. 2. Snapshots of surfaces obtained at (a) 6 = 0.5 (collapsed 
phase), (b) 6=1 (wormlike phase), (c) 6 = 30 (short linear 
phase), and (d) b= 120 (long linear phase). The surface size 
is A'" = 4842. The scale of the figures is difi'erent from to each 
other. 



First we show snapshots of surfaces of size N = 4842 in 
Figsl2l[a)-[2ljd) respectively obtained at 6 = 0.5 (collapsed 
phase), 6=1 (wormlike phase), 6 = 30 (short linear phase), 
and b = 120 (long linear phase). The scales are different 
to each other figure, because the size of surfaces extends 
from a small scale to a very large one as we will see in 
the mean square size later in this section. We find from 
Figsl2l[c) and[2ljd) that the surfaces in the short and the 
long linear phases are almost one-dimensional except the 
two terminal points, which have a two-dimensional struc- 
ture. 

Figures [3lja) and[3l^b) are the snapshot and the surface 
section obtained at &=140 (planar phase), and FigslSlJc) 
and[3]Jd) are those obtained at 6= 190 (spherical phase). 
The snapshots were drawn in the same scale, which is also 
identical to that of Figllja). From the section in FiglSl^b) 
of the planar surface, we understand that the surface is 
flat in the planar phase. It is also seen that the surfaces 




(c) 

Fig. 3. Snapshots of surfaces and their surface sections of size 
= 4842 obtained at (a),(b) b — 140 (planar phase) and at 
(c),(d) 6=190 (spherical phase). These figures were drawn in 
the same scale, which is also identical to that of Fig[2ja). 

are quite smooth in both phases. The smoothness is in 
sharp contrast to the case of the standard surface model 
such as the one in [27] , where the surface is not so smooth 
even in the smooth phase. 

The mean square size is defined by 



(3) 



where denotes the summation over the vertices i, and 
X is the center of the surface. The size and the shape of 
surfaces can be reflected in X^. 

Figures [4i;a),[4i;b), andH^c) show X^ vs. b of the 
1442 surface, the N = 2562 surface, and the N = 4842 
surface, respectively. Dashed lines drawn vertically in the 
figures denote the phase boundaries, where discontin- 
uously changes. The discontinuous changes of X^ indicate 
that the transitions are of first-order. The phase bound- 
ary between the collapsed phase and the wormlikc phase 
and that between the wormlike phase and the short linear 
phase are located in the small b region and therefore, they 
are not depicted in FigHl We can also see that X'^ wildly 
varies in the short linear phase on the surfaces of N — 2562 
and N = 4842 in FigslUb) and IH^c). X^ is sensitive to 
the size of the two-dimensional parts at the two terminal 
points of the surface in the short linear phase, and the size 
varies in that phase. For this reason, X'^ wildly varies in 
the short linear phase. 

In order to see the phase structure in the small b region, 
we plot in FiglS^a) X^ versus b obtained at this region. 
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Fig. 4. The mean square size versus b obtained on the surfaces of size (a) Ai' = 1442, (b) A'' = 2562, and (c) A'^ = 4842. Dashed 
hnes denote the phase boundaries. 
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Fig. 5. (a) The mean square size versus b in the small b region, (b) log-log plots of against A*' at the collapsed phase, 
at the planar phase and at the spherical phase, and (c) log-log plots of against TV at the wormlike phase and the long linear 
phase. 



The dashed lines were drawn at the phase boundary be- 
tween the collapsed phase and the wormlike phase and at 
the boundary between the wormlike phase and the short 
linear phase of the iV = 4842 surface. We find that con- 
tinuously varies at the boundary between the collapsed 
phase and the wormlike phase. This indicates that the 
transition between the two phases is not of first-order but 
of higher-order. On the contrary, the wormlike phase and 
the short linear phase are connected by a first-order tran- 
sition, because we see that discontinuously changes at 
the boundary. 

Figure E^b) shows log-log plots of X^ against N ob- 
tained in the collapsed phase at 6 = 0.5, in the planar 
phase at b close to the spherical phase, and in the worm- 
like phase at b close to the planar phase. Moreover, we 
show in Fig[5jc) log-log plots of X^ against N obtained in 
the long linear phase at b close to the planar phase, and 
in the wormlike phase at 6=1. 

The straight lines were drawn by fitting the data to 

X^ ^ N^/", (4) 

where H in the exponent is called the HausdorfF dimen- 
sion; the smoothness and the shape of surface can be re- 
flected in H. Thus, we have Hcoi (collapsed phase), H^or 
(wormlike phase), Hm (long linear phase), 77pia (planar 



phase), and Hgph (spherical phase), which are shown in 
Table [H 



Table 1. The Hausdorff dimension H obtained in the collapsed 
phase, in the wormlike phase, in the long linear phase, in the 
planar phase, and in the spherical phase. 



J/coi 11.2 ±0.89 H^or 1.85 ±0.11 
Hm 1.46 ±0.07 J/pia 1.89 ± 0.00 
gsph 1.99 ±0.01 



We note that iJgph = 1-99(1) is almost identical to the 
topological dimension of surfaces, and this result is consis- 
tent to the fact that the surface in the spherical phase is 
smooth and spherical as shown in FigElJc). On the other 
hand, we see that iJpia = 1.89(0) slightly deviates from the 
expected value H = 2, which is obtained under the uniform 
distribution of vertices on the disk. The reason of the de- 
viation of iJpia from H — 2 seems that the distribution of 
vertices is non-uniform as can be seen in the snapshot in 
FigEfa), where vertices are dense in the perimeter region 
while it is sparse in the central part of the surface. 
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H\\i — 1.46(7) also deviates from H = 1, which is the 
Hausdorff dimension of the straight line of uniform density 
of vertices. The surface in the long linear phase is not 
exactly identical to the straight line, because the surface 
includes the two-dimensional parts at the two terminal 
points. The surfaces in the wormlike phase are not exactly 
one-dimensional, and therefore the Hausdorff dimension is 
expected to deviate from H —I. This is a reason why the 
results i?wor = 1-85(11) deviates from the expectation. 

Finally, we comment on the Hausdorff dimension in 
the collapsed phase. The reason of Hco\ = 11.2(89) is larger 
than the physical bound H = 3 is that the surface is col- 
lapsed in the collapsed phase; the surface is allowed to 
self-intersect and is a phantom surface. If X'^ remained 
constant when N increased, H can be 7? — > cx). How- 
ever, in the model of this paper gradually increases 
with increasing N as we see in FiglS^b). This small and 
steady increase of X^ against N indicates that H has a 
finite value. We note also that the result i?coi > 3 is not 
always typical of the Nambu-Goto surface model. In fact, 
the Hausdorff dimension in the collapsed phase is larger 
than the physical bound; H > 3, in surface models with an 
intrinsic curvature energy [30] , where the surface tension 
energy is the Gaussian bond potential. 

The two-dimensional bending energy S3 is defined by 

S3^Y.^l~n,-n,), (5) 

where n.^ is a unit normal vector of the triangle i. We 
should note that 6*3 reflects the bending deformation, and 
therefore S3 changes its value depending on how wildly 
the surface fluctuates, although 5*3 is not included in the 
Hamiltonian. Therefore, the phase transition of shape trans- 
formation can also be characterized by anomalous behav- 
iors of 5*3, i.e., the discontinuous change of S3, because 
the shape transformation phenomena is expected to be ac- 
companied by the surface fluctuation phenomena as in the 
standard two-dimensional curvature model with Gaussian 
bond potential. 

Figure [6l[a)-[6jc) show the two-dimensional bending 
energy S3/NB versus b obtained on the surface of size 
N = 1U2, A^ = 2562, and A^ = 4842. Nb is the total number 
of bonds, which is given by Nb = 3N—6. The dashed lines 
were drawn in the figure at the same positions as those in 
Figs|Ha)-[4Kc) and therefore denote the phase boundaries. 
We find discontinuities of S3/NB at the phase boundaries 
in Figsinja), ini^b), andlHI^c). This observation also implies 
that the first-order transitions of surface fluctuation oc- 
cur at the phase boundaries indicated by discontinuous 
changes of in Figsl4l^a)-[4l[c). It is easy to understand 
that S3/NB — O in the spherical phase, because the surface 
is sufficiently smooth as shown in the snapshot in Fig[3fc). 
On the contrary, 1 — ■ rij has its maximal value at the 
edge of the disk in the planar phase, and 1 — • is 
expected be on the surface except at the edge. Thus, we 
understand the reason of the gap in S3 /Nb at the bound- 
ary between the planar phase and the spherical phase. 

The one-dimensional bending energy S2/NB versus b 
is shown in Figsl7^a)~(7l^c), where the surface size is = 



1442, iV = 2562, and A^ = 4842. The total number of com- 
binations 1 — t i • tj in S2 of Eq. ^ is identical to Nb the 
total number of bonds. For this reason, S2 is divided by 
Nb in Figsl7l^a)-l7l[c). The positions of the dashed lines are 
identical to those in the previous figures. The discontinu- 
ity is clear only at the boundary between the long linear 
phase and the planar phase. 

We also see that S2/NB continuously varies at the 
boundary between the planar phase and the spherical phase. 
This is inconsistent to the naive expectation that physi- 
cal quantities seem to have anomalous values at the edge, 
which is the singular part of the planar surface. The reason 
of the continuous behavior in 5*2 /Nb is because the ver- 
tices are dense in the perimeter region. Triangles become 
very thin at the edge of the planar surface, then 1 — • tj 
remains quite small at the edge and is comparable to that 
in the central part of the planar surface or in the spherical 
surface. 

Finally, we show the area energy Si/N in FigslSI^a)- 
[Hc). It is expected that Si/N = 3/2 from the scale in- 
variant property of the partition function. We find from 
Figsl5]Ja)-[5{c) that the relation Si/N = 3/2 is satisfied ex- 
cept in the long and the short linear phases. The surface in 
the linear phases is almost one-dimensional and therefore, 
the random shift of the vertices X in Monte Carlo updates 
seems not efficiently performed on such straight thin sur- 
faces. For this reason, Si/N seems to deviate slightly from 
3/2 in the linear phases. 

4 Summary and conclusions 

We have studied shape transformation phenomena in a 
surface model of Nambu and Goto, where the mechanical 
strength is given by a one-dimensional bending energy. 
The lattice is a triangulated mesh and is identical to the 
one widely used for surface models with two-dimensional 
bending energy. The Nambu-Goto surface model is known 
as an ill-defined model even with the standard two-dimensional 
bending energy. Thus, we aimed at showing that the model 
turns to be well defined when the one-dimensional bending 
energy is included in the Hamiltonian. 

It was shown that the model is well defined; there ap- 
pears no surface like a chestnut bur even at sufficiently 
small bending rigidity 6 — > 0. The surfaces have always 
a well-defined surface structure in the whole b region in- 
cluding 6— >0. 

We have also shown that the model has a variety of 
shapes including the smooth spherical surface and the col- 
lapsed surface. There are 6 different phases in the model, 
and almost all two neighboring phases are separated by 
discontinuous transitions. The appearance of the planar 
phase should be emphasized, because it is not yet re- 
ported that the fixed-connectivity surface model has a 
planar phase. It should also be emphasized that surface 
fluctuations are hardly seen in the inflated phases such as 
the planar phase and the spherical phase. The surfaces in 
those phases appear to be quite smooth. 

The results reported in this paper are very different 
from those in [35] , where the standard Gaussian bond po- 
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Fig. 6. The two-dimensional bending energy Ss/Nb versus b obtained on the surfaces of size (a) A'^ = 1442, (b) iV = 2562, and 
(c) A'^ = 4842. Dashed lines denote the phase boundaries. 
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Fig. 7. The one-dimensional bending energy S2/NB versus b obtained on the surfaces of size (a) AT = 1442, (b) iV = 2562, and 
(c) iV=4842. Dashed lines denote the phase boundaries. 
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Fig. 8. The area energy Si/N versus b obtained on the surfaces of size (a) iV = 1442, (b) iV = 2562, and (c) iV = 4842. Dashed 
lines denote the phase boundaries. 



tential is assumed as a bond potential. Bond potential 
considerably influences the phase structure of meshwork 
models. 
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